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^ . Abstract 

C ■ Using all the available empirical information, we analyse the spacing 

distributions of low-lying 2"*" levels in even-even nuclei by comparing them 
with a theoretical distribution characterized by a single parameter (the 
chaoticity parameter /). We use the method of Bayesian inference. We 
show that the necessary unfolding procedure generally leads to an over- 
, estimate of /. We find that / varies strongly with the ratio R4/2 of the 

excitation energies of the first 4"*" and 2^ levels and assumes particularly 
small values in nuclei that have one of the dynamical symmetries of the 
Interacting Boson Model. 

^ ! 1 Introduction 

The interplay between regular and chaotic motion in nuclei has been a 
k>( \ long-standing problem in Nuclear Physics. There is, on the one hand, 

j_j ■ overwhelming evidence in favour of simple dynamical models especially 

' in the ground-state domain. The evidence derives from the agreement 

between calculated and measured spectral properties. There is, on the 
other hand, equally strong evidence for the validity of a random-matrix 
description, especially from the spectral statistics of slow neutron reso- 
nances This success of random-matrix theory negates a dynamical 
description in terms of simple and (nearly) integrable models and has 
raised the question: Where in the spectrum of a nucleus with mass num- 
ber A does the chaotic region start? The statistical analysis of spectra 
needed to answer this question requires complete (few or no missing lev- 
els) and pure (few or no unknown spin-parities) level schemes. Some 15 
years ago, complete and pure level schemes were available for only a lim- 
ited number of nuclei (see, e.g., Refs. Elll])- The work of Ref. ^ then 
suggested that the nearest-neighbour spacing (NNS) distribution of low- 
lying nuclear levels lies between the Wigner and the Poisson distributions 
which are characteristic, respectively, of fully regular and fully chaotic 
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motion. Through the work of Refs. El El HI EOl Ell ED, the evidence 
presented in Ref. jS] has since become an established fact. 

The wealth of spectroscopic data now available in the Nuclear Data 
tables |13| has motivated us to investigate once again the nuclear ground- 
state domain. We are able to make more definitive and precise statements 
about regularity versus chaos in this domain than has been possible so far. 
As in Ref. we focus attention on 2+ states of select even-even nuclei. 
These nuclei are grouped into classes. The classes are defined in terms 
of the ratio J?4/2, i.e., the ratio of the excitation energies of the first 4"^ 
and the first 2^ level in each nucleus. We argue below that the classes 
define a grouping of nuclei that have common collective behaviour. The 
sequences of 2"*" states are unfolded and analysed with the help of Bayesian 
inference. The chaoticity parameter / defined below is determined for each 
class. The present paper summarizes two research papers |14l 115) where 
further details may be found. 

2 Data Set and Classification of Nuclei 

The data on low-lying 2"*" levels of even-even nuclei are taken from the 
compilation by Tilley et al. for mass numbers 16 < ^ < 20, from 
that of Endt |T7| for 20 < A < 44, and from the Nuclear Data Sheets E3 
for heavier nuclei. We considered nuclei for which the spin-parity as- 
signments of at least five consecutive 2"'"-levels are unambiguous. In cases 
where the spin-parity assignments were uncertain and where the most 
probable value appeared in brackets, we accepted this value. We termi- 
nated the sequence when we arrived at a level with unassigned , or when 
an ambiguous assignment involved a 2^ spin-parity among several possi- 
bilities, as e.g. = (2+, 4"*"). We made an exception when only one such 
level occurred and was followed by several unambiguously assigned levels 
containing at least two 2^ levels, provided that the ambiguous 2^ level 
is found in a similar position in the spectrum of a neighboring nucleus. 
However, this situation occurred for less than 5% of the levels considered. 
In this way, we obtained 1306 levels of spin-parity 2^ belonging to 169 
nuclei. The composition of this ensemble is as follows: 5 levels from each 
of 47 nuclei, 6 levels from each of 32 nuclei, 7 levels from each of 22 nuclei, 
8 levels from each of 22 nuclei, 9 levels from each of 16 nuclei, 10 levels 
from each of 14 nuclei, 11 levels from each of 5 nuclei, 12 levels from each 
of 2 nuclei, and sequences of 13, 14, 15, 17, 20, 21, 24, 30, and 32 levels, 
each belonging to a single nucleus. 

A class of nuclei is defined by choosing an interval within which the 
ratio 

i?4/2 = E{A+)/E{2+) (1) 

of excitation energies of the first 4^ and the first 2^ excited states, must 
lie. The width of the intervals was taken to be 0.1 when the total number 
of spacings falling into the corresponding class was about 100 or more. 
Otherwise, the width of the interval was increased. The use of the param- 
eter Q as an indicator of collective dynamics is justified both empirically 
and by theoretical arguments. We recall the reasons in turn. 
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(i) Casten et al. plotted -E(4^) versus E{2^) for all nuclei with 
38 < Z < 82 and with 2.05 < i?4/2 < 3.15. The authors found that 
the data fall on a straight line. This suggests that nuclei in this wide 
range of ^-values behave like anharmonic vibrators with nearly constant 
anharmonicity. As the ratio R4/2 approaches the rotor limit R4/2 ~ 3.33, 
the slope of the curve showing E{4:f) versus E{2f) decreases within a 
narrow range of i?(2^)-values, asymptotically merging the rotor line of 
slope 3.33. In a subsequent paper I19| it was found that a linear relation 
between E{4^) and E{2^) holds for pre-coUective nuclei with R4/2 < 2. 
Thus, from an empirical perspective, the dynamical structure of medium- 
weight and heavy nuclei can be quantified in terms of -R4/2. 

(ii) Theoretical calculations based on the Interacting Boson Model 
(the IBM-1 model \20\ ) support the conclusion that R4/2 is an appro- 
priate measure for collectivity in nuclei. The model has three dynamical 
symmetries, obtained by constructing the chains of subgroups of the U (6) 
group that end with the angular momentum group 50(3). The sym- 
metries are labeled by the first subgroup appearing in the chain which 
are !7(5), SU{3), and 0(6) corresponding, respectively, to vibrational, ro- 
tational and 7-unstable nuclei. Extensive numerical calculations for the 
classical as well as the quantum-mechanical IBM Hamiltonian by Alhassid 
et al. |21| indeed showed a considerable reduction of the standard mea- 
sures of chaoticity when the parameters of the IBM model approach one 
of the three cases just mentioned. The IBM calculation of energy levels 
yields values of R4/2 = 2.00, 3.33, and 2.50 for the dynamical symme- 
tries U{5), SU{3), and 0(6), respectively. Thus, we may expect increased 
regularity of nuclei having one of these values of i?4/2- 

One might expect that the chaoticity parameter / defined in Eq. ^ 
below also assumes small values for nuclei near magic numbers. For mass 
numbers in this domain, our data set is unfortunately too small to allow 
us to draw definitive conclusions. 

3 Statistical Analysis 
3.1 Chaoticity Parameter / 

To analyze the data, we need a guess for the form of the NNS distribution 
p{s, /). Here, s is the spacing of neighboring levels in units of the mean 
level spacing. The distribution p{s, /) depends on one or more parameters 
/ which describe the transition from Poissonian to Wigner-Dyson form. 
Several proposals have been advocated for p{s, /). Here we are guided by 
the following considerations. 

We consider a spectrum S containing levels which have the same spin 
and parity but may differ in other conserved quantum numbers which 
are either unknown or ignored. The A'-quantum number serves as an 
example. The spectrum S can then be broken down into m subspectra 
Sj of independent sequences of levels. Let fj,j = 1 . . .m with < fj < 
1 and fi ~ ^ denote the fractional level number, let Pj{s), j = 

1 . . . m denote the NNS distribution for the subspectrum Sj and p(s) the 
NNS distribution of S. Both p{s) and Pj{sj) are defined for spectra with 
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unit mean spacing. We assume that each of the distributions Pj{s) is 
determined by the Gaussian orthogonal ensemble (GOE). To an excellent 
approximation, the pj 's are then given by Wigner's surmise |23| 

Pw(s) = |sexp (-^s^) ■ (2) 

The construction of p{s, f) for the superposition is due to Rosenzweig 
and Porter |24| . It depends on the (m — 1) unknown parameters fj, 
j = 1, (m — 1). This fact poses a difficulty because in practice, we do 
not know the composition of the spectrum. We are not even sure of how 
many quantum numbers other than spin and parity are conserved. To 
overcome the difficulty, we use an approximate scheme first proposed in 
Ref. |25|. Effectively, we replace the (m — 1) parameters fj by a single 
one, the mean fractional level number / = fj. This leads to an 
approximate NNS distribution for S, 

P{sj) = [l - / + / (0.7 + 0.3/) ^ 

xexp|-(l-/)s-/(0.7 + 0.3/)^| . (3) 

We use / as a fit parameter. 

For a large number m of subspectra, / is of the order of 1/m. In this 
limit, p{s,f) approaches the Poisson distribution as it should. On the 
other hand, when / ^ 1 the spectrum approaches the GOE behaviour as 
it must. This is why we refer to / as to the chaoticity parameter. If the 
spectrum S is not pure but rather a superposition of subsequences corre- 
sponding to different values of an ignored or unknown quantum number 
then the mean value / of the fractional density of the superimposed se- 
quences is smaller than unity, and the composite sequence looks rather 
like a sequence of levels with mixed dynamics. 



3.2 Unfolding 

Prior to the actual statistical analysis, every sequence of levels has to be 
unfolded |22| to obtain a new sequence with unit mean level spacing. In 
the case of a single long spectrum, unfolding is a standard procedure. It 
consists in fitting a slowly varying function e{E, a) to the experimental 
staircase function N{E) of the integrated level density. The fit is obtained 
by optimizing a set of parameters a. The function e depends monotoni- 
cally on the energy E. Therefore, we can transform E to e. With respect 
to the new energy variable e, the level density is uniform and equal to 
unity. 

If the available ensemble of spacings consists of many short sequences 
of levels (we call this a "composite ensemble" ) , unfolding is not standard 
nor is it altogether irrelevant. To test the standard unfolding procedure, 
we have generated short sequences of levels from three artificial ensembles 
containing 50, 100, and 200 spacings. Construction of the latter involves 
an artificially chosen chaoticity parameter /o and is described in the fol- 
lowing paragraph. These are referred to as the "initial" ensembles. Each 
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short sequence is then artificially folded with a monotonically increasing 
function of energy. An unfolding procedure is subsequently applied to 
each sequence. The unfolding procedure does not trivially reproduce the 
initial ensembles and yields the "final" ensembles. The chaoticity param- 
eter / is then determined for the final ensembles using a x'^ fit and the 
Bayesian method described below. 

The ensembles of spacings are constructed with the help of a random- 
number generator. We choose average spacing unity and /o = 0.6 for the 
chaoticity parameter. This value is close to what has been obtained in 
the previous analysis |S] of low-lying nuclear levels. We generate a set of 
spacings that obeys the probability distribution with f ~ fo- In this 
way, we generate three "initial" artificial ensembles of 50, 100, and 200 
spacings. Our procedure is open to the criticism that our construction 
does not pay attention to the stiffness of GOE spectra. We are in the 
process of rectifying this shortcoming. 

The test of the unfolding procedure leads to the following conclusions, 
(i) Using several unfolding functions leads to nearly the same values for 
/. This confirms the insensitivity of the final ensemble of spacings to the 
form of the unfolding function, (ii) The unfolding procedure introduces a 
bias towards the GOE, i.e. the best-fit value of / is larger than fo- This 
is borne out by both, the Bayesian inference and the x^-analysis of the 
spacing histograms for the final distributions. The trend increases as the 
lengths of the short sequences is decreased. This is simply understood: 
The unfolding of sequences of just two levels each would give a delta- 
function peaked at the value of unity (the mean level spacing) and, thus, 
show strong preference for the GOE. The trend becomes weaker as the 
sequences become longer but disappears only in the limit of very long 
sequences. As a consequence, the analysis of the nuclear data set will 
reliably yield only relative values of /. 

The actual unfolding of the data was done by fitting a theoretical 
expression to the number N{E) of levels below excitation energy E. The 
expression used here is the constant-temperature formula 



The three parameters A'o, Eq and T obtained for each nucleus vary con- 
siderably with mass number. Nevertheless, all three show a clear ten- 
dency to decrease with increasing mass number. For the effective temper- 
ature, for example, we find, assuming a power-law dependence, the result 
T = (15 ± 4)^-(o e2±o o5) rpjjjg ^g^jyg jg consistent with an analysis 

of the level density of nuclei in the same range of excitation energy carried 
out by von Egidy et al. |4j. These authors find T = (19 ± 2)>l-(0 '38±0 02) 
MeV. 

3.3 Bayesian Analysis 

Given Eq. ^ for the proposed distribution, we apply Bayesian analy- 
sis to the data. Let s = (si, S2, sjv) denote a set of spacings Sj. We 
take the experimental spacings Sj to be statistically independent. This 
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assumption does not apply in general. Indeed, the GOE produces signif- 
icant correlations between subsequent spacings. However, we recall that 
we are interested only in the NNS distribution. This distribution is only 
weakly affected by correlations. We calculate the posterior distribution 
for / given the events s. We first determine the conditional probability 
distribution p (s |/) of the set of spacings s = (si, S2, sat) for a fixed /. 
We accordingly write 

JV 

P{^\f)-l[p{s.j) , (5) 

i = l 

with p{si, f) given by Eq. Bayes' theorem then provides the posterior 
distribution 

of the parameter / given the events s. Here, fi{f) is the prior distribution 
and ^ 

M(s)= / p{s\f)^^if)df (7) 



is the normalization. We use Jeffreys' rule |2(i| 



Kf) 



p{s\f) [d\np{s\f)/dfY ds 



1/2 

2 



(8) 



to find the prior distribution. The latter can be interpreted as the distri- 
bution ascribed to / in the absence of any observed s. It is approximated 
by 

fi{f) = 1.975-10.07/+48.96/^-135.6f^-f205.6/''-158.6/^+48.63/ . 

(9) 

Even for only moderately large N, it is useful to write p{s\f) in the form 

p(s|/) = e-'^^'^' , (10) 

where 

= (l-,f)(s> + ^/(0.7+0.3/)(s^)-(ln[l-/+J/(0.7+0.3/)s]) . (11) 

Here the notation (x) = (1/A'') ^^^^ a;^ has been used. By calculating 
the mean values (• ■ •) in Eq. lllll for various spectra, one finds that the 
function <^(/) has a deep minimum, say ed f — f^. One can therefore 
represent the numerical results in analytical form by parametrizing as 

^if) = A + B{f - hf + C{f - hf . (12) 

We then obtain 

P(/|s) = c^^if) eM-N[B{f - fif + C{f - hf]) , (13) 
where c = e-^^/M{s) is a normalization constant. 
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The last step ol the Bayesian analysis consists in determining the best- 
fit value of the chaoticity parameter / and its error for each NNS distri- 
bution. When P (/ |s) is not Gaussian, the best-fit value of / cannot be 
taken as the most probable value. Rather we take the best-fit value to 
be the mean value / and measure the error by the standard deviation aof 
the posterior distribution 0, i.e. 



The results obtained for / and a are given in Figure 1 of Ref. [H]. Fig- 
ure 2 of that reference shows a comparison of the spacing distributions 
conditioned by / and the histograms for each class of nuclei. In view 
of the small number of spacings within each class, the agreement seems 
satisfactory. 

We recall that the analysis of many short sequences of levels tends to 
overestimate /. Therefore, we focus attention not on the absolute values 
of / but on the way / changes with -R4/2. The graph of / against R4/2 in 
the Ref. |14| has deep minima at R4/2 = 2.0, 2.5, and 3.3. These values 
of R4/2 are associated with the dynamical symmetries of the Interacting 
Boson Model mentioned above. Another minimum of statistical signifi- 
cance occurs for 2.25 < R4/2 < 2.35. This minimum may indicate that 
nuclei which lie between the limiting cases of the (7(5) and 0(6) dynami- 
cal symmetries, are relatively regular. One may associate this region with 
the critical point of the ?7(5)-0(6) shape transition in nuclei. lachello )27l 
has recently shown that this transition is approximately governed by the 
"critical" E{5) dynamical symmetry. Nuclei with E{5) dynamical sym- 
metry have R4/2 ~ 2.2. Experimental examples of this critical symmetry 
have been found by Casten and Zamfir |28l . 

In summary, we have determined the chaoticity parameter / for 2^ 
levels of even-even nuclei at low excitation energy with the help of a sys- 
tematic analysis of the NNS distributions. While in a single nucleus the 
number of states with reliable spin-parity assignments is not sufficient 
for a meaningful statistical analysis, a combination of sequences of levels 
taken from similar nuclei provides a sufficiently large ensemble. As the 
measure of similarity we have taken the ratio R4/2 of the excitation ener- 
gies of the lowest 4^ and 2^ levels in each nucleus. The mean chaoticity 
parameter / is found to be indeed dependent on -R4/2. It has deep minima 
at R4/2 — 2.0, 2.5, and 3.3. These minima correspond, respectively, to 
the ?7(5), SO{Q), and 5'C/(3) dynamical symmetries of the IBM. A further 
minimum may relate to the critical £'(5) symmetry. 




(14) 



This is not optimal but provides a useful approximation. 



4 Results and Discussion 
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